Energy dissipation statistics in a shell model of turbulence 
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The Reynolds number dependence of the statistics of energy dissipation is investigated in a shell 
model of fully developed turbulence. The results are in agreement with a model which accounts 
for fluctuations of the dissipative scale with the intensity of energy dissipation. It is shown that 
the assumption of a fixed dissipative scale leads to a different scaling with Reynolds which is not 
compatible with numerical results. 
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t-H . One of the most important problems in fully developed turbulence is the description of the energy transfer mech- 
anism. In stationary situations, the energy injected at the large scale £o, transfers at rate e down to the dissipative 
CD ■ scale 77, where it is removed at the same rate by viscous dissipation. The fundamental assumption in the study of 
fully developed turbulence is that in the limit of very high Reynolds numbers Re, the energy dissipation e becomes 
independent of Re (i.e. of the viscosity, being Re = u^l^fv, with uq a typical large scale velocity) Jl] ||. In the same 
limit, the Kolmogorov theory predicts universal scaling of the velocity structure functions in the inertial range of 
scales rj < £ < £o- 

s q (e) = ((Su(£)) q )^ui (£) <q (1) 

o : 

Q\ with exponents ( q = q/3. 

Several decades of experimental and numerical investigation have shown that scaling laws (Q) are indeed observed 
but with exponents £ g corrected with respect to the Kolmogorov prediction || . This is the essence of the intermittency 
problem, which has received a lot of attention in the modern approach of the study of fully developed turbulence. 
Experiments have shown that intermittency also affects energy dissipation statistics || which is not uniform in 
i ^T *' the turbulent domain. A phenomenological description of intermittency is the multifractal model [ ^] . This model 
1 ■ introduces a continuous set of scaling exponents h which relate the velocity fluctuations entering in (n ) with the large 
O \ scale velocity u : 

Su(£)~u [-) . (2) 

; f e \z(h) 

The exponent h is realized with a probability \ j^J where Z(h) is the codimension of the fractal set on which the 
/i-scaling holds. The scaling exponents of structure functions (0) are obtained by a steepest descent argument over 



exponents h: 



( q = inf [ph + Z{h)\ . (3) 

h 



The scaling region is bounded from below by the Kolmogorov dissipation scale rj at which dissipation starts to 
dominate, i.e. the local Reynolds number is of order 1: 

At variance with the Kolmogorov theory, in the multifractal description of intermittent turbulence, the dissipative 
scale is a fluctuating quantity. This implies a series of consequences which have been investigated in past years [^,^) . 
As we shall see later the description of the fluctuations of the dissipative scale is crucial for the correct evaluation of 
the Reynolds number dependence. 

In this Paper we are interested in the dependence of the statistics of energy dissipation on the Reynolds number. 
The physical picture is that dissipation becomes more and more intermittent as the Reynolds number increases. 
Assuming that the multifractal description can be pushed down to the dissipative scale, one predicts for the moment 
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of energy dissipation a power-law dependence on Re, with exponents related to the structure function exponents ([j]) 
pO|,p[. We will see that this prediction is rather natural and confirmed by numerical simulations on a shell model. 

The dimensional argument for the prediction goes as follows. In a dimensional approach, the energy dissipation is 
evaluated as 

From (^) and (Q) one has that 77 ~ igRe^ 1 ^ . Inserting in (||) and computing the average of the different moments, 
one ends with the expression 

dfi(h) Re - e p Re~ 0p (6) 

where the integral has been evaluated by a steepest descent argument (assuming Re — > 00) and 

~3ph-p + Z(h)~ 
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The standard inequality in the multifractal model (following from the exact result £3 = 1), Z(h) > 1 — 3h, implies 
for (0) 9(1) = which is nothing but the request of finite nonvanishing dissipation in the limit Re — > 00. For p > 1, 
p < 0, i.e. the tail of the distribution of e becomes wider with Reynolds number. 

Let us stress that the above argument is only a reasonable dimensional argument. It is essentially based on two 
assumptions: a physical one concerning the fluctuations of the dissipative scale according to (^), and a more formal 
one on the possibility of extending the multifractal description down to the dissipative scales. The two assumption 
are independent: indeed, as we will see, it is possible to give different predictions by changing assumption (ji|) |Tl| . 

It would thus be important to address the problem with experiments or direct numerical simulation at high Reynolds 
numbers. Recent high resolution DNS gives support to (]?]) Q, but the Reynolds number is not large enough to 
discriminate clearly between different predictions. 

Shell models are extremely simplified models of turbulence. Nevertheless, they are deterministic, nonlinear dynam- 
ical models which display intermittency and anomalous scaling exponents reminiscent of real turbulence [ p"3[ . Their 
main advantage is that with shell models one can perform very accurate simulations at very high Reynolds numbers; 
for this reason they are thus natural candidates for a numerical investigation of Reynolds number dependence. 

In shell models, the velocity fluctuations are represented by a single complex variable u n on shells of geometrically 
spaced wavenumber k n — ko 2". The particular model we adopt for our investigation is a recently introduced model 
which displays strong intermittency corrections fl4| . The model equations are 

= ik n \ U n+ 2U* n+x - ^U n+1 U* n _ 1 + ~li„_lU„_2^ - vk^u n + /„ (8) 

where v is the viscosity and /„ is a forcing term restricted to the first two shells. For v = f n = the model 
conserves the total energy E = J2 n \ u n\ 2 - For simplicity, the forcing adopted for the present simulations is f n oc , 
which guarantees a constant energy input e. The large scale Reynolds number of the simulation is estimated as 
Re — e 1 / 3 / (vk^ 3 ) and is numerically controlled by the value of the viscosity. 

The chaotic dynamics is responsible for intermittency corrections to the structure functions exponent C?j here 
defined by means of 

S 9 (n) = (|u„H ~ k~ c « , (9) 

which are close to the experimental values |fl3| . In Figure [l] we plot the spectrum of structure function exponents 
obtained from very long simulations. The multifractal codimension Z(h) is numerically obtained from £(p) by inverting 
the Legendre transform (^|). The result is shown in Figure p|. We observe that, because of the strong intermittency in 
the model, it is numerically difficult to obtain statistical convergence of structure functions (||) of order q > 8. As a 
consequence, the minimum exponent for Z(h) is /i mm ~ 0.2. 

From the energy balance equation we have the instantaneous energy dissipation 



(10) 



2 



which average is (e) = e in stationary conditions. 

We have performed very long simulations at different Reynolds numbers, starting from Re = 2 x 10 5 up to Re = 
10 8 . For each simulation we computed the different moments of energy dissipation, (e p ). Shell models dynamics is 
characterized by strong bursts of energy dissipation which limits the possibility of computing with confidence high 
order moments. Here we limited to moments p < 8. In Figure || we plot the behavior of (e p ) as a function of Re 
for different values of p. The power law behavior is evident and the scaling exponent 9 p can be estimated with good 
accuracy. By construction (e) = e is independent of Re. 

The scaling exponent 9 P (JsJ) are plotted in Figure |], together with the multifractal prediction (Q) . Let us observe 
that, because the largest q in (||) is q — 8, the estimate of p in (0) is limited to values less than p ~ 2.5. For higher p, 
the numerical evaluation of (Q) feels the effect of the cutoff of h on h m i n . Nevertheless, we have a rather large range 
of moments (0 < p < 2.5) over which numerical data display a perfect agreement with (Q). 

As discussed above, prediction (||) makes use of the fluctuating dissipative scale r\. If, on the contrary, one assumes 
that dissipation scale enters in (||) as an averaged quantity the prediction for 9 P is different: assuming v (du(fj) 2 ) / fj 2 ~ e 
as the definition of the (non-fluctuating) dissipation scale fj (this is the only choice that ensures that (e) ~ Re ) , one 
ends up with 9 p = [p( 2 - £^]/[2 - £,] fl- 

Our results allow us to discriminate between the prediction (^) and the one obtained with a non-fluctuating dissi- 
pative scale. Figure ^ shows that the numerical 9 P is definitely not compatible with the latter alternative, whereas it 
supports with good accuracy the prediction (Q). 

In conclusion, it is an expected consequence of the existence of intermittency in the energy transfer that the 
dissipative scale r\ fluctuates according to the local intensity of energy dissipation, being smaller where the dissipation 
is stronger and vice versa. The fluctuations of the inner scale of turbulence reflect onto the Reynolds dependence 
of the statistics of energy dissipation. Long numerical simulations of shell-models confirm with great accuracy the 
validity of the multifractal model, which accounts for the fluctuations of rj, and rule out alternative models which do 
not describe properly the correlations between r] and e. 
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FIG. 1. Structure function exponents £ g obtained from simulation of the Shell Model. The number of shells is N = 24 and 
v — 10~ 7 corresponding to Re = 10 s . 

FIG. 2. Z{h) computed by inverting the Legendre transformation from the data of Figure |l]. 

FIG. 3. Moments of energy dissipation (e p ) as function of Re for p — 1 (+), p = 2 (x) and p — 3 (*). 

FIG. 4. Energy dissipation scaling exponents P . Symbols represent the exponents obtained from the fit of Figure []. 
Continuous line is prediction (Q) taking into account the fluctuations of the dissipative scale. Dashed line is the prediction 
obtained by assuming an average dissipative scale. 
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